Influence of Quantum Fluctuations on Phase Coherent Andreev Tunnehng 
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We study the subgap transport properties of a small capacitance normal metal-superconductor 
tunnel junction coupled to an external electromagnetic environment. Mesoscopic interference be- 
tween the electrons in the normal metal strongly enhances the subgap conductance with decreasing 
bias voltage. On the other hand, quantum fluctuations of the environment destroy electronic phase 
coherence and suppress the subgap conductance at low bias (Coulomb blockade). The competition 
between charging effects and mesoscopic interference leads to a non-monotonic dependence of the 
differential subgap conductance on the applied bias voltage. This feature is pronounced, even if the 
coupling to the environment is weak and the charging energy is small. 
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Charge transport through a tunnel barrier between a 
normal metal (N) and a superconductor (S) is a widely 
investigated topic ^ . At energies much smaller than the 
superconducting gap A, tunneling of single particles is 
exponentially suppressed. Under these conditions, charge 
transport through an N-S interface is dominated by An- 
dreev reflection Q. If the normal metal and the super- 
conductor are separated by a low transparency tunnel 
barrier, two-electron tunneling [|j determines the subgap 
conductance. Recently, there has been much interest in 
the subgap properties of mesoscopic N-S junctions 
The dependence of the subgap conductance of a tunnel 
junction on temperature or applied bias voltage can be 
very different, depending on the precise lay-out of the 
system under consideration. 

If, for instance, we consider a tunnel junction between 
a superconductor and a thin metallic film at low tem- 
peratures, electrons move phase coherently in N and 
undergo multiple elastic scattering events by impurities 
or rough sample boundaries. As a consequence, they 
will be scattered back to the junction interface several 
times, where they attempt to tunnel into S. Two-electron 
tunneling involves two almost time reversed electrons. 
Therefore, the phase of the two-electron tunneling am- 
plitude is not randomized by elastic scattering and the 
amplitudes for various tunneling attempts add up coher- 
ently. This strongly enhances the subgap conductance at 
low bias voltages [^,0]. On the other hand, in mesoscopic 
N-S tunnel junctions with a small junction capacitance 
C, charging effects Q become important. In order to 
tunnel, the two electrons should overcome the charac- 
teristic Coulomb interaction energy Ec = e^/C. This 
will strongly suppress the subgap conductance [^JlO| at 
low temperatures and bias voltages fcsT, eV <C -Ec, a 
phenomenon known as Coulomb blockade of two-electron 
tunneling. 

In the present paper we will discuss the influence of the 
competition between charging and interference effects on 
the subgap conductance of a single N-S tunnel junction. 



Charging effects in a single junction are conveniently de- 
scribed using the so-called electromagnetic environment 
model In this model, electron tunneling is studied 
in the presence of quantum phase fluctuations due to the 
Johnson-Nyquist noise of the external circuit, seen by the 
junction. In the simplest case, this circuit consists of a 
capacitor C {i.e., the junction capacitance) and an exter- 
nal series resistor Z{u!) = R, see Fig. pi. The influence 
of such an environment on the subgap conductance of 
N-S junctions has been studied before |l2|, but without 
considering mesoscopic interference. We will show how 
interference effects are destroyed as the series resistance 
R, which determines the coupling to the environment, 
is increased. If the charging energy is smaller than the 
superconducting gap, the competition between quantum 
fluctuations and mesoscopic interference leads to a non- 
monotonic dependence of the differential subgap conduc- 
tance on bias voltage, even for weak coupling. 

The system shown in Fig. |l|a can be described by the 
Hamiltonian H = + Ht. Here, Hq denotes the un- 
perturbed Hamiltonian, Hq = Hn + Hg + H^nv , where 
Hiq and Hs describe the disordered normal metal and 
the superconductor, respectively. The electromagnetic 
environment is described by the usual bosonic Caldeira- 
Leggett [ p^ Hamiltonian TJenw ■ The tunnel Hamiltonian 
Ht transfers electrons between N and S and couples the 
electrons to the environment; it will be treated perturba- 
tively. In the interaction picture Ht takes the form 
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Here, ^i^cr is a fermionic field operator for an electron 
with i = N, S and spin a =|, J,; T(r, r') is the amplitude 
to tunnel from a point r in N to r' in S, and V the applied 
bias voltage. The phase operator ip{t) describes the volt- 
age fluctuations at the tunnel junction, induced by the 
electromagnetic environment. Its dynamics is governed 

by Henv 
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From an expansion in Ht up to fourth order, using lowing expression for the subgap current: 
standard imaginary-time techniques, one obtains the fol- 
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xCAr(ri,r'2,r3,r4;0,Ti,T2,r3)$(0, ri,T2,T3)Cs(r4,r3,r'i,r'2;T3,r2,ri,0) , 



(2) 



where /3 = l/ksT and the analytical continuation for the 
bosonic Matsubara frequencies iwi — > —eV/h+i6, iuo2 
eV/h + iS, iuj^ — > eV/h + i6 has to be performed. In 
order to obtain (Q), we assumed tunneling to occur be- 
tween neighboring points on the barrier B, located at 
zb — 0; correspondingly we put T{r,r') = tod^{r — 
r')d{z). The amplitude tg can be expressed in terms 
of the normal state conductance Gt of the tunnel bar- 
rier, Gt = ATrtlSBNp/ihvFRK)- Here Rk = 2nh/e^ 
is the quantum resistance, Np the density of states at 
the Fermi level and vp the Fermi velocity; Sb is the 
area of the barrier surface. Furthermore, we introduced 
the four-point correlator (r i , r2 , r3 , r4 ; Ti , r2 , r3 , r4 ) = 
(f^*|^(ri,Ti)*|^(r2,r2)*|,(r3,r3)*^,(r4,T4)) describ- 
ing the propagation of two electrons for i = N, S, as 
well as a four-point phase correlator $(0, ri, T2, r3) = 
(^f^^~^Ma)+v{Tl)^'P{r2)-p{T3)]) related to the voltage fluc- 
tuations. The averages (. . .) are taken with respect to 
i/AT, Hs, and Henv, respectively; in addition the correla- 
tor C is to be averaged over disorder. 

Further simplification can be achieved following Ref. 
1^. At energies much smaller than the gap A, two elec- 
trons propagate coherently through N over distances of 
the order = y/ hD / m.ax[eV, kBT, h/Tip], much larger 
than the corresponding length ^5 = y/hD/A in S {D is 
the diffusion constant and the phase breaking time); 
moreover, the lifetime ~ h/A of a quasiparticle in S in 
the intermediate state is negligibly small. Therefore, we 
have Cs ~ S{ri - r2)S{r3 - ri)S{Ti)SiT2 - T3) in Eq. (|). 
The dominant contribution to the subgap current can 



now be written as 
2>RkG^ 
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The integrand of Eq. (||) is depicted diagrammati- 
cally in Fig. |^. We see two electrons tunneling from 
S to N at initial position and time (0;0) thereby in- 
teracting with the environment. The electrons propa- 
gate through N to position r, where they arrive at time 
t; their propagation is described by the disorder av- 
eraged two-particle correlator (Coopcron), CN^r^r) — 
(CAr(0, 0, r, r; 0, 0, r, r))disordcr. Then they tunnel back 
into S, interacting once more with the environment. The 
wavy line in Fig. |^b denotes the phase correlator 4'(t) — 
(TVe~2«[v(r)-v(o)]^ = exp[4J(T)], where we choose J ac- 
cording to the electromagnetic environment model 
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Here Zt{iS) = l/(icjC+ is the total impedance 

seen by the junction. 

Upon performing the analytic continuation in Eq. (|^) 
the pair tunneling current finally is found to be 



1{V) = 



3Rk G^ 
2TT^eSBNp Jb 



f d^r f dEdE' CN{r]E)P(E') 



1 - exp[-2ey/3] 
1 - exp[(£;' - 2eV)(3] 



ifHE - 2eV + E')/2] - f[{E + 2eV - E')/2]). 



(5) 



The function CN{r; E) is the spatial fourier transform of 
the real part of the diffusion propagator l/[—iE+fiDQ^+ 
h/Tip]. The probability P{E) to emit or absorb a photon 
with frequency E/h during tunneling is defined as 



PiE)^^ J ^^exp [AJ{t) 

where J{t) can be obtained from Eq. 
for r > 0. 



-iEt/h], 



(6) 



by putting t = it 



As an example, we will study the Andreev current 
of a N-S tunnel junction, consisting of a quasi one- 
dimensional (ID) normal metal wire in contact with a su- 
perconductor via a tunnel barrier with dimensions much 
smaller then ^jy. For a quasi ID wire with cross section 
5vK, we have 
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in Eq. (|5|). The junction has a capacitance C and is em- 
bedded in a purely resistive environment Z{uj) = R. The 
total impedance is thus given by Zt(w) = 1 / [iuoC + R^^) . 
However, the calculation of P{E) can only be per- 
formed numerically in this case jl4|. Therefore, in order 
to proceed analytically, we will use the approximation 
Re[Zt((jj)] = _Rexp(— w/wc), which has the correct zero- 
frequency limit Re[Zt(0)] = R. The cut-off frequency 
is chosen to be lOc — EcRk /^R, such that the approx- 
imated impedance also gives the correct phase correla- 
tion function J{t) ~ exp{—i2Ect/h) in the limit of short 
times. In particular, this guarantees that this approxima- 
tion yields the correct behavior for i? = oo, namely the 
complete suppression of tunneling for eV < Ec = e^/C. 

At zero temperature the phase correlation function is 
readily found to be <&(<:) = l/ll + iujct]"; its fourier trans- 
form (||) yields the probability distribution 



P{E) 



-E/huJc / ^ 
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The negative part of the spectrum is truncated, since 
at T = photons can only be emitted. The power 
a = 8R/Rk can be interpreted as a parameter determin- 
ing the coupling strength between the electronic phase 
and the environment. 

For the non-interacting system, a = 0, the subgap con- 
ductance is proportional to the "coherence resistance" 
Rcoh = (.n/o'Sw Hi, where a denotes the conductivity 
of the wire. In the limit of strong coupling, a = oo, a 
gap appears in the TV curve below Ec- We will focus 
on the case of small charging energies Ec A. For a 
finite coupling and finite phase coherence time [ p^ , 
two regimes exist (see Fig. ||): (i) For very small voltages 
eV <C h/r^p, interference is cut off by r^p and the I-V curve 
shows a power law behavior / cx . This is what one 

would expect for noninterfering electrons, (ii) For higher 
voltages eV ^ h/r^ the coherence length is voltage de- 
pendent and the I-V characteristic changes: / cx l/^+^Z^. 
Remarkably, if = oo and a = ac = 1/2, suppression 
of the current by charging effects and enhancement by 
interference cancel each other exactly at voltages below 
Ec/e, such that the TV curve is linear. For values of a 
larger than the "critical" coupling the power is al- 
ways larger than one and a Coulomb gap starts to evolve 
with increasing a. 

The differential conductance G{V) = dl/dV is 
strongly suppressed at eV < h/T^ for arbitrary a > 0. 
The zero bias peak, which is the fingerprint of phase co- 
herent Andreev tunneling, is destroyed by charging ef- 
fects. Instead, the differential conductance will display a 
peak at finite bias (Fig. |^) . Increasing the voltage on the 
one hand lifts the Coulomb blockade, but on the other 
hand decreases the coherence length. If a < Uc, the max- 
imum in the differential conductance appears at a bias 
eV ^ h/r^ and will shift to zero bias for = oo. The 



coupling to the quantum fluctuations of the environment 
is too weak to fully destroy phase coherence. Therefore 
coherent pair tunneling is blocked only at voltages be- 
low h/r^, where interference is cut off. For a > ac, the 
coupling is strong enough for charging effects to domi- 
nate the behavior. The conductance is suppressed in the 
entire voltage regime below Ec by quantum fluctuations, 
which makes Tip superfluous as a cutoff for the divergence 
of G at zero bias. The maximum appears at eV Ec 
and shifts towards Ec as a is increased. 

At finite temperatures, the electrons can gain en- 
ergy by absorbing environmental modes. Therefore, the 
Coulomb blockade is gradually lifted with increasing tem- 
perature. The probability distribution P{E) at small 
energies \E\ <C htOc can be calculated analytically, fol- 
lowing p^ ]. In the long time limit, the phase correlator 
reads $(i) ~ [(TrfesT/nwc)/ sinh(7rtfcBr/n)]". From (|) 
we find 



P{E) 



exp(£'/2fcBT) 
2T:T{a)uJc 



2-KkBT 



huJr 



iE 
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For low temperatures ksT <^ hujc one easily calculates 
the zero bias differential conductance G{T, V = 0) with 
the help of (|). For fc^T < ?i/t^ we find G(T, F = 0) - 
T", whereas G(T, V = 0) ^ ya-i/2 foj. ^^j^ ^ ^/^^^ 

In order to determine G{T, V) away from zero bias, 
the function P{E) should be calculated for arbitrary E. 
For nonzero temperatures, this can only be done numer- 
ically. Qualitatively, one expects the Coulomb blockade 
to be lifted if ksT > h/T,p for weak coupling, a < 
whereas for a > ac thermal smearing becomes relevant 
at temperatures fcsT ^ Ec- As an example, the differen- 
tial conductance as a function of bias voltage is sketched 
in Fig. ^ for a < ac at various temperatures . 
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FIG. 1. (a) Single N-S junction with a capacitance C cou- 
pled to an external circuit with an impedance Z(lo) and a 
voltage source, (b) Two electrons propagating coherently in 
a disordered normal metal as a cooperon C(r;r) (half-moon) 
coupled to the electromagnetic environment by the phase cor- 
relator ^(t) (wavy line). The upper loop describes two elec- 
trons which immediately form a Cooper pair after entering 
the superconductor; the lower loop describes the correspond- 
ing time-reversed process. 
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FIG. 2. Andreev current in units of Jo = {3Gj'/2TvaSwe) 
X ^JhDEc for T = and ft/r^ = 0.5-Ec. Curves from top to 
bottom correspond to a = (dashed line), ct = 1/4, 1/2, 2, 8 
(solid lines) and a = oo (dashed line). 
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FIG. 3. Diff erential c onductance in units of Go = 
{'iGT/'2%aSw)yJ{'hD/Ec) for T = and h/r^ = 0.5Ec. The 
maximum evolves to the right as a is increased from a = 
(dashed line), taking a = 1/8, 1/4, 1/2, 2, 8 (solid lines). 



1 




FIG. 4. Sketch of differential conductance for a < ac- 
From top to bottom, temperature increases: T — (dotted 
line), fcsT ^ h/r^p (solid lines), and ksT < h/r^p (dashed 
lines) . 
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